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; Abstract 

The exclusive decay of T to a vector plus pseudoscalar charmonium is studied in perturbative 
^ QCD. The corresponding branching ratios are predicted to be of order 10 for first three T reso- 

nances, and we expect these decay modes should be discovered in the prospective high-luminosity 
e + e" facilities such as super B experiment. As a manifestation of the short-distance loop contri- 
bution, the relative phases among strong, electromagnetic and radiative decay amplitudes can be 
deduced. It is particularly interesting to find that the relative phase between strong and electro- 
magnetic amplitudes is nearly orthogonal. The resonance-continuum interference effect for double 
r~| . charmonium production near various T resonances in e + e~ annihilation is addressed. 
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I. INTRODUCTION 



Very rich J/ip decay phenomena have historically served as an invaluable laboratory to 
enrich our understanding toward the interplay between perturbative and nonperturbative 
QCD [l|, 0]. By contrast, much fewer decay channels of T are known to date. It would be 
definitely desirable if more knowledge about bottomonium decay can be gleaned. 

The typical branching fraction for a given hadronic decay mode of T is in general much 
smaller than that of J/if>. It is partly due to the smaller QCD coupling at the b mass scale 
than at the c scale, and more importantly, it is because the branching ratio gets diluted by 
a scaling factor of (m c /rrib) n when descending from charmonium to bottomonium (here n is 
some number no less than 4). These might intuitively explain why very few exclusive decay 
modes of bottomonia have been seen so far. 

Due to rather large b mass, T not only can dematerialize into light hadrons, it also can 
decay to charmful final states. In this work, we plan to study the exclusive decay of T 
into double charmonium, or more specifically, J/ip(ip') plus r] c (r]' c ). The hard scales set by b 
and c masses in this type of processes justify the use of perturbative QCD (pQCD). Since 
the involved mesons are all heavy quarkonium, it is natural to employ NRQCD factorization 
approach [3]. This work constitutes a continuation ofprevious studies on bottomonium decay 
to double charmonium, namely, Xb, Vb ~ ¥ J/^J/^ Although these decay modes have 

not yet been seen, some experimental information have already been available for inclusive 
J / ip production rate from T decay (1, 0, S] : 

B[T{1S) -> J/ip + X] = (6.5 ±0.7) x 1(T 4 , 

B[T(lS)^i/>' + X] = (2.7 ±0.9) x 1(T 4 , 

B[T(25) -> J/i/j + X] < 6 x 10- 3 , 

B[T(4S) -> J/ip + X] < 1.9 x 10" 4 . (1) 

These inclusive decay ratios set upper bounds for our exclusive processes. It is worth noting 
that T — > J/iprjc violates the hadron helicity conservation 0, T(|. It is thus natural to expect 



that the corresponding branching fractions are very suppressed. 

One important impetus of this work is from the double charmonium production at T(AS) 



resonance measured by Belle in 2002 11]. The observed cross section is usually entirely as- 
cribable to the continuum contribution because of rather broad T(4S*) width. Nevertheless 
for a full understanding, it is worth knowing precisely the impact of the resonant decay on 
the measured double charmonium cross section. Furthermore, stimulated by Belle's discov- 
ery, a natural question may arise-what is the discovery potential for double charmonium 
production in e + e~ experiments operated at lower T peaks? Since the first three T reso- 
nances are much narrower than T(4S), the resonant decay contribution should dominate 
over the continuum one. Our study is motivated to answer this question. 

One interesting problem in exclusive decays of a vector quarkonium is to know the relative 
phase between strong and electromagnetic amplitudes. For example, the corresponding 
relative phase in J/ip — > PV (P, V stand for light h and 1 mesons) has been extensively 
studied and found to be nearly orthogonal [Hf, [H, 14, [H, 17]. In our case, the relative 
phase naturally emerges as a short- distance effect and thus is perturbatively calculable. 
Curiously, it is also found to be approximately orthogonal. 

The rest of the paper is organized as follows. In Section [Til we present the lowest-order 
NRQCD calculation for the decay process T — > J/ip + r) c , including strong, electromagnetic 
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and radiative decay channels. In Section IHIl we present the predictions to the branch- 
ing fractions for various T decays to double charmonium, and conclude that the discovery 
potential of these decay modes is promising in the prospective Super B experiment. We 
also discuss the relative phases among three amplitudes, putting particular emphasis on the 
nearly orthogonal relative phase between strong and electromagnetic amplitudes. The con- 
nection between our results and the previous discussions on the nearly 90° relative phase in 
J/ip decays is remarked. In addition, we also study the impact of the resonance-continuum 
interference on J/ip + 7] c production cross sections at various T resonances in e + e~ exper- 
iments. We summarize and give a brief outlook in Section IIVI In the Appendixes, we 
illustrate how to analytically derive some loop integrals that appear in Section [Til 



II. COLOR-SINGLET MODEL CALCULATION 

The process T — > J/ip + r] c can proceed via three stages: the bb pair first annihilates into 
three gluons, or two gluons plus a photon, or a single photon; in the second step, these highly 
virtual gluons/photon then convert into two cc pairs, which finally materialize into two fast- 
moving .S-wave charmonium states. Due to the heavy charm and even much heavier bottom, 
both the annihilation of bb and creations of cc pairs take place in rather short distances, it 
is thereby appropriate to utilize pQCD to study this hard exclusive process. 

This process is somewhat similar to the widely studied J/tp — > PV decay, but bears 
the virtue that applicability of pQCD should be more reliable. It is commonly believed 
that some nonperturbative mechanisms should play a dominate role in many charmonium 
exclusive decay processes, where the credence of pQCD seems rather questionable. This 
consensus is exemplified by the notorious pir puzzle [2J, |l8 . 

While it is customary to use the light-cone approach to deal with hard exclusive processes 
involving light mesons (for a recent attempt to study J/ip — > pir from this perspective, see 
Ref. [19j), it is for our purpose most proper to employ an approach embodying the non- 
relativistic nature of quarkonium. NRQCD factorization is a widely accepted effective-field- 
theory framework to describe the quarkonium inclusive production and decay processes, 
which incorporates systematically the small velocity expansion j3j. Although a rigorous 
formulation for exclusive quarkonium decay has not yet been fully achieved within this 
scheme, one may still be well motivated to work with models akin to the NRQCD ansatz. 

The color-singlet model can be viewed as a truncated version of NRQCD approach, in 
which one still assumes a factorization formula, i.e., the decay rate can be separated into 
the perturbatively calculable part and universal nonperturbative factors, however only with 
the contribution from the color-single channel retained. We do not know how to include the 
possible color-octet contributions in a clear-cut way, but it is plausible to assume their effects 
are unimportant for reactions involving only S'-wave quarkonium as in our case. Notice that 
NRQCD and color-singlet model are often referring to the same tool in literature, so we will 
also use them interchangeably. 

Let Q, P and P signify the momenta of T, J/ip and t] c , respectively. In color-singlet model 
calculation, one starts with the parton process b(pb)b(pi) — > c(p c )c(pc) +c(p c )c(p E ), then 
projects this matrix element onto the corresponding color-singlet quarkonium Fock states. 
This work is intended only for the zeroth order in relativistic expansion, hence we can neglect 
the relative momenta inside each quarkonium, i.e., set Pb = Pi = Q/2, p c = Pc = -P/2, and 
Pc = Pc = P/2. For the bb pair to be in a spin-triplet and color-singlet state, one simply 
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FIG. 1: Some representative lowest-order diagrams that contribute to T — > 3g — > J/ip + ?? c - 



replaces the product of the Dirac and color spinors for b and b by the projection operator 



1 



2rrih 



? T X 



' ^t(O) J ® 1 



For the outgoing J /if} and r] c , one makes the following replacements: 



(2) 



v(Pc) u{Pc) 

v(p B ) u{p c ) 



' fj^(P + 2m c ) x ( -1-^(0) ) ® 1 



2^2 
1 

2^ 



i7 5 (f + 2m c ) x 



1 \ 1 



(3) 
(4) 



Here Ey and £j,^ are polarization vectors for T and J/if>. N c = 3, and l c stands for the unit 
color matrix. The nonperturbative factors t/>t(0), ipj/^(0) and Vv(0) are Schrodinger wave 
functions at the origin for T, J /if) and rj c , which can be inferred either from phenomenological 
potential models or extracted from experiments. By writing (J2J), (EJ) and (HI) the way as they 
are, it is understood that M T = 2m b and Mj/^ « M^ c = 2m c have been assumed. 

Before moving into the concrete calculation, we recall first that since both strong and 
electromagnetic interactions conserve parity, the decay amplitude is then constrained to have 
the following Lorentz structure: 



M = Ae^E^E* J %Q a P /3 . 



(5) 



Apparently, J /if; must be transversely polarized in T rest frame. All the dynamics is encoded 
in the coefficient A, which we call reduced amplitude. Our task in the remaining section 
then is to dig out its explicit form. 



A. Three-gluon Amplitude 

We begin with the strong decay amplitude. Some typical lowest-order diagrams are shown 
in Fig. [H which starts already at one loop order. Using the projection operators in (j2J), ([3]), 
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and dH), we can write down the corresponding amplitude: 



M 3s = 2 t,(7-7*7«) t r(T°{T\ T*}) ^ ( °> ^Mg> J_ J_ 



x 



tr[(ff + 2m b ) gVyT(fe + m 6 ) 7 "] tr[(@ + 2m b ) fal v {- fa + m b )^Y] 

^2 — m b ^1 ~ m 6 

m fe tr[$ + 2m b ) fa + m b )Y(fa + m b )Y] 



(kj - m 2 h ) {k\ - m 2 b ) 

tr^* (P + 2m c ) 1 ^ c - fa + m c ) lvlB {f + 2m c ) lp \ 
x " 7Z u m —2 > ( 6 ) 



(p c - &i) 2 - m 



where two internal gluons carry momenta k\ and &2, respectively, which are subject to 
the constraint fa + k 2 = ®. Some elaboration is in order. Because T has charge conju- 
gation quantum number —1, three intermediate gluons must arrange to the color-singlet 
state <i a6c |a)|6)|c) [d abc (f abc ) represents the totally (anti) symmetric structure constants of 
SU(N C ) group]. This restriction removes all the possible 0(g®) diagrams involving 3-gluon 
vertex. As a result, we only need retain those Abelian diagrams in which each of three 
gluons is connected between the b and c quark lines in both ends. There are totally twelve 
such diagrams, but it turns out that for each of diagrams, there is another one generating 
exactly identical amplitude, which explains the prefator 2 in the right hand side of ([6]). 
Among the six diagrams needed to be considered, one can further divide them into two 
groups: one carries a color factor oc ti(T a T b T c )tr(T a T c T b ) , whereas the other carries that 
oc tr (T a T b T c )tr (T a T b T c ) . These two groups yield identical reduced amplitudes except this 
difference. Thus we only need consider three diagrams with distinct topologies, as depicted 
in Fig. [TJ and incorporate the following color factor: 



, abc _ (iV c 2 -l)(iV c 2 -4) 



tr(T a T b T c ) tr(T a {T b , T c }) = - d abc d abc = ^ ^ >- , (7) 



which reassures us that only those intermediate gluons with overall C = — 1 can contribute 
to this process. 

Straightforward power counting reveals that the loop integrals in ([6]) are simultaneously 
ultraviolet and infrared finite. In absence of the need for regularization, we have directly 
put the spacetime dimension to four. 

After completing the Dirac trace in (jUJ), we end up with terms in which the Levi-Civita 
tensor is entangled with the loop momentum variable. Since all these terms will finally 
conspire to arrive at the desired Lorentz structure as dictated in (jSJ), we may exploit this 
knowledge to get rid of the antisymmetric tensor prior to performing the loop integral [20] . 
First we may identify the partial amplitude M^ v through M. = M pu e*f^. Equation (j5J) 
then demands 

= Ae^ aP Q a P^, (8) 
Contracting both sides of (jSj) with e^ upa Q p P a , one can extract the reduced amplitude using 
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where |P| = [(Q ■ P) 2 — Q 2 P 2 ] 1 / 2 /My is the modulus of the momentum of J/ip (r] c ) in the 
T rest frame. 

After this manipulation is done, we end in a concise expression 



A 



3g 



2^2(N 2 C -1)(N 2 C -A) nc4 



5/2 



m 



7/2 I 



Vt(0)Vj/v»(0)^ c (0)/ 



m: 



mt 



where / = f x + f 2 + / 3 , and 



./i 

/3 



(i 4 ^ (m 2 - 4m 2 ) (A; 2 - m 2 ) + h ■ (3Q - P) h ■ P - (1 + m 2 c /ml)(h ■ Qf 



Z7T 



k 2 k 2 (k 2 - ml) (k 2 -h-P) 



d 4 h (m 2 - 4m 2 ) (k 2 -ml) + k 2 ■ P k 2 ■ P - (m 2 /m 2 ) (k 2 ■ Q)< 



m 2 h 



k\ ik\ - m 2 b ) k\ (k\ -h-P) 
d A h h-(Q- 2P)(k 2 -h-P)- 2(m 2 - 4m 2 )fc! ■ k 2 



k\ (k 2 - m 2 b ) k 2 (k 2 - ml) (k\ -h-P) 



(10) 

,(H) 
(12) 

(13) 



Since /j is dimensionless, it can depend upon m^ and m c only through their dimensionless 
ratio ml /ml. These loop integrals can be worked out analytically, and the results are 
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(14) 



(15) 



where Li 2 is the dilogarithm function, and (3 = y/1 — 4£. We will illustrate in Appendix |A] 
how to obtain this result. The emergence of imaginative part of / characterizes the contri- 
bution from two on-shell internal gluons. The shapes of the real and imaginary parts of / 
are displayed in Fig. [31 

It is instructive to know the asymptotic behavior of / in the £ — > limit. This can be 
readily read out from (1T41) and (1T51) . 
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FIG. 2: Three representative lowest-order diagrams that contribute to T — > gg'j — > J/ip + rj c , where 
the J /if) comes from the photon fragmentation. 



Note both the real and imaginary parts blow up logarithmically in the limit £ — > 0, as 
can be clearly visualized in Fig. [31 These (quadratically) logarithmical divergences in the 
m c — > limit are obviously of infrared origin. Nevertheless, this does not pose any practical 
problem, since a nonrelativistic description for a zero-mass bound state, as well as the 
resulting predictions, should not be trusted anyway. It is interesting to note that, provided 
that £ is not overly small, say, £ > 10 -4 , then — Im/ is always bigger than |Re/|, or more 
precisely phrased, — < arg/ < — | . 



B. Two-gluon-one-photon Amplitude 

We next turn to the contribution through the radiative decay channel. C-parity conserva- 
tion demands that one end of photon line must be attached to the b quark. Those diagrams 
obtained from replacing one gluon by one photon in Fig. [1] do contribute, however their 
magnitudes are much less important than the diagrams shown in Fig. [2, which essentially 
proceed as T — > gg(-+ r) c ) +7(— * J ftp)- This is because in the latter case, the J /if} is created 
via the photon fragmentation, which thereby receives a ml/ ml enhancement relative to the 
former. We will only consider the latter case, in which the lowest-order contribution also 
starts at one loop. Using the projection operators in (j2J), (jSJ), and (@J, it is straightforward 
to write down the corresponding amplitude: 

Vt(0)Vj/^(0)^(0) fd% 1 1 



ka _ o AT-1/2 fr(T a T b^ fr(T a T b^ p p 2 4 ^(0)^(0) %M f 

M gai - 2 N c tr(T T ) tr(T T ) e b e c e g s — J ^ ^ fc 



4 U2 L2 



X 



+ 



+ 



tr[(# + 2m b ) fa fj/^P- A + m b )Y(fa- fa + m b )^] 
((P -p b ) 2 - mg)((p 6 - h) 2 - ml) 

tr[($ + 2m b ) faYih~ A + rribh^A- P + m b ) fj^\ 
((P-p b ) 2 -m 2 )((p b -k 2 ) 2 -m 2 ) 

tr[(# + 2m b ) faYifa- A + m b ) fj/^b- fa + m b )Y 
({Pb - ki) 2 - m 2 )((p b - k 2 ) 2 - ml) 



x tr[7 5 (f + 2m c )7 /J (^ c - fa + m c ) lv \ 
(p c - h) 2 - m 2 

where the momenta carried by two internal gluons are labelled by ki, k 2 , which satisfy k\ + 
k 2 = P. The factor 2 in the right side of (fT8l) takes into account the identical contributions 
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FIG. 3: Real and imaginary parts of /(£) and <?(£). 



from other three crossed diagrams. 

Following the same shortcut adopted in the 3g channel, we can derive the desired reduced 
amplitude with recourse to Eq. (jUJ), 



&y/2(N? — 1) e h e c ^aa 2 mV 2 , , . , , . , , . f m 2 \ 

a„ = ^ J (m g t : 2m / ) | P ; 2 ,„g ^(0)^,(0)^(0) s {-£ ) , (w) 

where the dimensionless function g is defined by 

/ m 2 c \ rd 4 ki (2m 2 c /mlQ ■ r - P ■ r) P ■ r + 2(ml - 4m 2 c )r 2 
9[ — )= j— k%k%(kl - h ■ Q)(k% - k 2 ■ Q) ' m) 



mi I .1 itc 



For convenience, we have introduced a new internal momentum variable r, which is defined 
through h\ = P/2+r and k 2 = P/2 — r. Note that the integrand is symmetric under r — > — r, 
reflecting the symmetry k\ <-> k%. A gratifying fact is that the charm propagator has now 
been cancelled in the denominator. We dedicate Appendix [B] to a detailed derivation of 
this loop integral. Like its counterpart / in the three-gluon channel, the function g is both 
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ultraviolet and infrared finite. Its analytic expression reads 



(1 - 20 ln[2 - 4£] + 4v^(Rj tan" 1 J ^ 
e/3 |4tanh- 1 /3 ln[2£] + 2Li 2 [-/3] - 2Li 2 [/3] + Li 2 
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l + /3) 2 + 4zv^(r^O 



Li 2 
Li 2 



(l-/3) 2 + 4zy^(r^O 

4(1 - 20 
"/3(l+/3) 2 + 4z/3v^(r^O 



Li 9 



4(1-20 

/3(l-/3) 2 + 4z/y0T^Q 
4(1 - 20 



4(1 - 20 
27r£/3tanlT 1 /?. 



(21) 
(22) 



The shapes of real and imaginary parts of g are displayed in Fig. [3] Note that — Reg is 
always bigger than — Img for any £, or put in another way, — 7r < argg < — Apparently, 
the imaginary part of g vanishes as £ — > 0, whereas the real part of g approaches the following 
asymptotic value: 
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Re<?(0 = - T + ln2 + O(£ln0 



(23) 



In contrast to /, both of the real and imaginary parts of g admit a finite value in the £ — > 
limit. Not surprisingly, the asymptotic behavior of this function is quite similar to the 
analogous one in the T — > rjcj process (20| . 



C. Single-photon Amplitude 



Let us now consider the electromagnetic contribution via the annihilation of bb into a 
single photon, with some typical diagrams shown in Fig. HI This process is closely related to 
the contiuum J/ip + r] c production in e + e~ annihilation, which has recently aroused much 
attention since the measurements were first released by Belle collaboration [ll|. Rather 
unexpectedly, it is shortly found that the leading-order NRQCD prediction to the pro- 
duction cross section falls short of the data by about one order of magnitude [2l|, 22], 
which su bsequen tly trig gered intensive theoretical efforts to resolve this alarming discrep- 
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ancy 

In the Born order, one can directly import the time-like electromagnetic form factor of 
S-w&ve charmonium first deduced in Ref. 2l| to here, and the corresponding lowest-order 
one-photon amplitude reads 



16v / 2(A r2 - 1) 7r 2 efee c a a s 
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FIG. 4: Two representative lowest-order diagrams that contribute to T — ► 7* 
are totally four diagrams in class a) and two in class b). 



J l 'ijj + rj c . There 



where the second term in the parenthesis represents the pure QED contribution in which 
J /if) arises from photon fragmentation, as is represented by Fig. HJd). 

Recent calculations indicate that the J/i[) + r\ c electromagnetic form factor is subject to 
large perturbative and relativistic corrections at B factory energy [28, 30]. It seems that 



the disturbing discrepancy between B factories measurements and NRQCD predictions have 
been largely resolved once these large corrections are taken into account. Motivated by this, 
from now on we will replace the entities in the parenthesis in (1241) by a positive constant K 
(> 1), which presumedly encompasses all the radiative and relativistic corrections. 



D. Decay Width and Asymptotic Scaling Behavior 

It is now the time to lump three different contributions together. Plugging (TTUj) . (JUJ), 
and (T2"4~]) into the formula 

|p|3 

T[T -> J/i) + 77J = i-i- \Ay + A 3g + A ggi \ 2 , (25) 

lZ 7T 

we then obtain the desired decay partial width. Note the cubic power of momentum reflects 
that Jfij) and n c are in relative P wave. This formula has already taken into account the 
spin average of T and the polarization sum over J/ip. The result is 

r[T^ J/tP + vd = r[T^e+e-] 



2 20 7r 2 e 2 a 2 IPI 3 
x Q Tji 



T 



where a 7 = K, 



«3 9 = -=^-^/(4), (27) 
72iiebe c a |P| 2 \ m 6/ 



a 991 — A „ /_2 0™2M1DI2™2 9 [ o ' ( 2 ^) 



m 2 



4-7T (m 2 — 2m 2 )|P| 2 m 2 
and r[T — >■ e + e _ ] = lQirela 2 ^^) / M\ is the electronic width of T. 
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It is instructive to deduce the asymptotic behaviors of these three different contributions. 
Because we are more concerned about the power-law scaling, we will take f,g ~ 0(1) for 
simplicity since they vary with quark masses logarithmically at most. Assuming ^jm(O) ~ 
~ ( m c^c) 3 ^ 2 {v c is the typical relative velocity between c and c), from (1261) we find 



r[T 




-> W + »7c] 




r[T - 


* e+e~] 


r[T 




-> W + 




r[T - 


+ e+e~] 


r[T - 




-> J/V> + ?7 C ] 




r[T -» 


e+e~] 



~ (30) 

or m£ 

~«^*e 8 - (31) 



First interesting observation is that both (129!) and (1301) exhibit the \jm\ scaling behavior. 
This is as expected from the celebrated helicity selection rule inperturabative QCD, which 
is applicable for both single-photon and three-gluon processes [9|. The reason is as follows. 
The final-state J/if> must be transversely polarized, in line with the parity and Lorentz 
invariance, the hadron helicity conservation Aj/^ + X Vc = is violated by one unit, hence the 
ratio is suppressed by an extra 1/Vn| relative to the leading-twist 1/mf scaling. In contrast, 
the corresponding ratio in ggj channel, (13T1) . though suppressed by coupling constants with 
respect to other two subprocesses, nevertheless enjoys a much milder (~ kinematical 
suppression, because the J/ ip directly comes from the photon fragmentation. Simple power 
counting implies that these three different contributions have comparable strengths for the 
physical masses of b and c. 

Another noteworthy fact is that, there are relative phases among three amplitudes, which 
are encoded in the / and g functions. Since all these phases originate from the loop integrals, 
we may regard them of short-distance origin. 

III. PHENOMENOLOGY 

A. Determination of K from B factories measurement 

First we want to determine the value of K in (1261) . which characterizes the magnitude 
of higher-order corrections to the single-photon amplitude. For the sake of simplicity, we 
will assume the K factors are equal in our case and in J /if} + r] c production through e + e~ 
annihilation to a virtual photon. Of course, this is just an approximation, because the virtual 
gluon line connecting b quark and final-state c quark, as well as the relativistic correction in 
T, which will emerge in our process accounting for the radiative and relative corrections, are 
absent in the double charmonium production in continuum. We will assume these additional 
corrections are insignificant. 

First we recall the continuum double charmonium cross section in the lowest order in a s 
and v 2 c 0,0: 



2 2 vr 2 2 2 |p| 3 

<w[e + e- - J/i, + Vc ] = a, +tl - 1#/*(0)<(0) , (32) 

where 07^- = ^f-- For simplicity, the pure QED contribution where J/ip is produced via 
photon fragmentation (the analogous diagram to Fig. Hb) has been neglected. 
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TABLE I: Experimental inputs for T(nS) and S-wave charmonia (taken from Ref. 3l|)- The last 
column lists the wave functions at the origin for various S-wave charmonium states, retrieved from 
the measured electric width through f|33|) by assuming m c = 1.5 GeV and a s (2m c ) = 0.26. 



u 


Mass ^Vjev J 


r tot (kev J 


J-e+e- (.KeVJ 


ipH\y) (LreV ' J 


T(15) 


9.460 


54.02 ± 1.25 


1.340 ±0.018 




T(25) 


10.023 


31.98 ±2.63 


0.612 ±0.011 




T(35) 


10.355 


20.32 ± 1.85 


0.443 ± 0.008 




T(45) 


10.579 


20500 ± 2500 


0.272 ± 0.029 




J/1> 


3.097 




5.55 ±0.14 


0.263 


Vc 


2.980 






0.263 




3.686 




2.48 ±0.06 


0.176 


rfc 


3.638 






0.176 



In this work, we extract the wave functions at the origin for vector charmonium states 
from their measured electric widths. We will use the formula incorporating the first order 
perturbative correction 

rf ,/, + -i 4:ire 2 c a 2 2 / 8a a (2m c ) \ 2 

T[J/ij -> e+e ] = -—^-^(O) ^1 — ] • (33) 

Heavy quark spin symmetry is then invoked to infer the wave functions at origin for the 
corresponding 1 Sq charmonium states. All the involved charmonium wave functions at origin 
are tabulated in Table [H 

If we choose m c = 1.5 GeV, a s = 0.22, we then obtain from (1321) the tree level continuum 
J/ijj + r] c cross section at y/s = 10.58 GeV to be 4.74 fb. This theoretical prediction can be 



contrasted with the most recent B factories measurements [32j, |33j: 

a{e + e~ -> J/tp + r} c ] x B^ 2 = 25.6 ± 2.8(stat) ± 3.4(syst) fb , Belle 
a[e + e- -> J/ip + rj c ] x B% = 17.6 ± 2.8(stat)^;f(syst) fb , BABAR (34) 

where B^ 2 represents the branching ratio of r) c decay to more than 2 charged tracks, hence 
should be less than 1. With large uncertainties, both measurements seem to be marginally 
consistent with each other. 

If we assume the measured o" cont [e + e~ — > J/ip + r) c ] to be 23 fb, and expect that large 
radiative and relativistic corrections to (|32p can bring the leading-order NRQCD prediction 
to this value, we then require K = a/23/4.74 w 2.2. This K factor is r oughly consistent with 
what is obtained through actual higher-order NRQCD calculations [28j, |30j. Although we 
extract this constant through the T(AS) — > J ftp + rj c process, we will assume it is universal 
in all other double charmonium decay channels of T(nS). 



B. Exclusive decay of T(nS) to double S- wave charmonium 

To date, T exclusive decays to double charmonium have not yet been experimentally 
established. To make concrete predictions from (126]) . we need specify the values of all the 
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TABLE II: Predicted partial widths and branching ratios for various decay channels of Y(nS) to 
vector plus pseudoscalar charmonium. 



Decay channels T (eV) }3 


Decay channels T (eV) 13 


T(15) -> J /if) + Vc 0.208i°;^2 3.9+6-6 x 10 -6 


T(2S) J/^ + r? c 0.08218:5^8 2.6±f;g x lO" 6 


T(lS)-> J/V^ + ^c 0.109±g;J|I 2.0±?"| x 10~ 6 


T(2S)-J/^ + t/ c 0.042l8;8E? 1.3ljj : * x 10" 8 


T( 15)^ V ,/ + r?c 0.0931°;^ 1.7l^xl0- 6 


T(25) V' + % 0.037l8:8i 1-lla? x 10 ~ 6 


T(15)-^' + V c 0.045i°;°30 0.8±J;f x lO" 6 


t(25)^v' + ^ c o.oi7i8:8 2 ? o.5i8ixio- 6 


T(35) -> J/V> + Vc 0.054±g;g|| 2.7l?;g x 1(T 6 


T(45) - J/V + % 0.03118:8^ 1.5±g;| X 10~ 9 


T(3S)-J/^ + t/ c 0.02718;^ 1.3±g;J x 10" 6 


T(4S)^J/^ + V ' c 0.01518^ 0.7l5| xio- 9 


r(3S)^V' + % 0.024l°-8?4 1.2l^xl0- 6 


T(45)^V' + r? c 0.01418:85s o.7i5 : ° x io- 9 


T(35)-^ + i/ c O-Olllg-E? 0.6181 xlO- 6 


t(4S)-v" + »/ c 0.00718:851 o.3l8^xio- 9 



input parameters. We fix m c to be 1.5 GeV, but take m& as a variable- for each T(nS) 
decay process, we approximate it as half of My( n s)- The magnitude of |P| is determined 
by physical kinematics. We assume K = 2.2 for all decay channels, and take the values of 
the wave functions at the origin for various charmonium from Table [B As for the coupling 
constants, we take a = 1/137, and cY s (m&) = 0.22. The uncertainties of our predictions are 
estimated by sliding the renormalization scale from 1m\, to m^jl (corresponding to varying 
a s from 0.18 to 0.26). It should be cautioned that the ambiguity of the inputted b mass, 
especially for higher T excitations, can bring even more severe uncertainty due to the higher 
powers of m& appearing in (12"6"1) . 

Our predictions to the partial widths and branching ratios for all decay channels are listed 
in Table [TTJ One clearly sees that the branching fractions for all decay processes (except 
for T(4S")) are about 10 _ 6 ; which are perfectly compatible with the measured inclusive J /if) 
production rates from T(nS) decay, Eq. p]). It is interesting to note that our hadronic 
decay processes have even smaller branching ratios than the radiative decay T — > 7^7 {B ~ 
3 x IO" 5 ) |20|. This may be partly understood by that T[T -> r] cl }/T[T -> e+e~] ~ 

which has a milder 1 / ml scaling behavior compared to the 1 / suppression in our processes, 
as manifested in Eqs. (I29I) and (I30I) . 

Between 2000 and 2003, CLEOIII has recorded about 20 million, 10 million and 5 million 
decays of T(IS'), T(2S) and T(3S), respectively [sj]. So there should be a few to tens of 
produced events for each double charmonium mode. Unfortunately, because the cleanest 
way of tagging J / if) is through the dimuon mode, only 6% fraction of the produced events 
can be reconstructed. Further taking into account the acceptance and efficiency to detect 
/i, it seems rather difficult to observe these double charmonium production events based on 
the existing CLEOIII data sample. By contrast, the high luminosity e + e~ colliders such as 
Belle and BABAR have already collected a enormous amount of data at T(AS) peak. If 
they could dedicate some significant period of run at the lower T resonances, it is feasible 
for them to discover these decay channels unambiguously. Needless to say, the discovery 
potential is very promising for the planned super-high-luminosity e + e~ facility like Super B 
factory. 

It is important to understand the interference pattern among three different amplitudes. 
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In our case, the phase in each amplitude manifests itself as short- distance effect arising from 
loop, and is perturbatively calculable. Let us take T(1S) — > J/if> + r] c as an example. Taking 
f = 4m2/M| « 0.10 and a s = 0.22, we find from ([27]) and (125} 



3.89 e 



-il05° 



0.44e l24 °. 



(35) 



Curiously, the strong decay amplitude is almost orthogonal to the electromagnetic amplitude, 
while the radiative decay amplitude is almost in phase with the electromagnetic one. It is 
also obvious to see that the strong decay amplitude has the most prominent strength, the 
electromagnetic one the next, and the radiative decay amplitude the least. 

In digression, it may be instructive to know the relative strengths of three different 
channels in inclusive T decay. From the following experimental inputs: 



R 



Ri, 



R~ 



rpr -/*+/*-] 

r[T - ggg] 



r[T — > 

r[T -» ggj\ 
^^ggg] 



= 3.56 ±0.07, 
39.11 ±0.4, 



34 



35[, 



0.027 ±0.003, 



36 



(36) 



we can infer 



B[T -> ggg] : B[T 



7 



X] : B[T -> ggi\ = 82.7% : 7.5% : 2.2% , 



(37) 



where these three branching ratios sum up to 1 — B[T — > = 92.5%, as they should 1 . 
A very simple expectation is that each amplitude in an exclusive process scales with the 
corresponding v^mci- The relative strengths of three amplitudes in fl35l) roughly respect 
this scaling rule if one assumes K — 1. Nevertheless, the truly important point is that, the 
orders of strengths of three amplitudes are same for both inclusive and exclusive decays. 

We can gain more intuition about the interference pattern by examining the individual 
contribution to the partial width. Had we retained only a 7 in ( |26i) . the partial width for 
T(IS') — > J/ij) + T] c would be only 0.065 eV. If we kept a 3g only, the width would instead be 
0.204 eV. If we include both a 7 and a^ g but discard a ggi , the width would become 0.210 eV, 
which is rather close to the full answer listed in Table HTl 0.208 eV. This numerical exercise 
clearly corroborates our expectation about the relative importance of these three different 
channels. 

The phase structures in (|35|) also hold for other decay channels of T(nS) to double 
charmonium. We take T(4S") — > J/ip + r] c as second example to verify this point. Taking 
£ = Am 2 c /My^ s) 0.08, we obtain 



a 3g 



4.20 e 



-i 102° 



"991 



0.52e 



i20° 



(38) 



It has been of great interest to decipher the interference pattern between the strong and 
electromagnetic amplitude in J /if) decays. The relative phase between 3q and 7 amplitude 
in J /if) -> PV has been determined to be around -(106 ± 10)° [H, M, Q H EE E3- This 



1 We have not included the contribution from the radiative transition T 
negligible branching ratio. 



?7b7, which has a completely 



14 



is surprisingly close to our rinding in T decay. Suzuki has argued that the large relative 
phase in J/ip decay must arise from long-distance rescattering effect, and emphasized that 
it is impossible for the perturbative quark-gluon process to generate it [16j . However, our 
calculation provides an explicit counterexample against his claim, showing that the short- 
distance contribution alone suffices to generating such a large relative phase. 

It is worth mentioning that some years ago, Gerard and Weyers argued there should be 
universal orthogonality between strong and electromagnetic amplitude for each J/ip exclu- 
sive decay mode [37[ • This assertion may seem to be backed by numerous phenomenological 
evidences 2 . They have attributed this orthogonality simply to the orthogonality of gluonic 
and one photon states. Inspecting their arguments carefully, one finds that they only prove 
the incoherence between three-gluon and single-photon decays at inclusive level, whose va- 
lidity crucially relies on summing over all possible decay channels. Since there is no room 
for such a summation for exclusive J/ ip decay, there is no any simple reason to believe why 
strong decay amplitude should be orthogonal to the electromagnetic amplitude channel by 
channel. 

Because their reasoning is based on rather general ground, one may test it in T exclusive 
decay. As a matter of fact, we can directly present a counterexample. Imagine a fictitious 
world with an extremely heavy b quark, say m;, ~ Mpi an k, but with an ordinary charm quark. 
For the would-be T decay to J /if) + r) c , we then find from (HM and (|T7|) that the phase of / 
is very close to zero, so is the relative phase between a^ g and a 7 . 

One may wonder why Gerard and Weyers's assertion seems to enjoy considerable success 
when applied to J/ip decays, even though it looks theoretically ungrounded. One possible 
explanation is that, due to some specific dynamics, the relative strength and phase between 
electromagnetic and strong amplitudes are roughly identical for each J/ip exclusive decay 
mode, preserving the same pattern as in the inclusive decay. This approximate scaling 
between exclusive and inclusive channels is exemplified in the discussion following fl37j) . 
This pattern does not necessarily hold for other vector quarkonium decays. 

It is straightforward to see that, the approximate —90° phase between strong and elec- 
tromagnetic amplitude in our process is simply a consequence of the not-too-tiny mass ratio 
m c/ m 6 ~ 0.1 and the opposite electric charges of c and b (see left panel of Fig. [3] and (j27j) ). 
It may seem to be a marvellous coincidence that the relative phase determined in our case 
is very close to that in J/ip — > PV, especially regarding that the latter process should be 
largely dictated by nonperturbative long-distance dynamics. We don't know exactly which 
nonperturbative mechanism should be responsible for the universal orthogonal phase in var- 
ious J/ip decay modes. It is fun to notice that, however, in the constituent quark model, the 
masses of u, d and s quarks are several hundreds of MeV, consequently d s /m 2 c rs m^/m 2 , 
so our formalism seems to be able to explain the nearly orthogonal phase in J/ip — > PV 
entirely within the short- distance quark-gluon picture. 

Lastly we stress that the phases determined in (15BT) and (J2HD are subject to large un- 
certainties. Since they are determined only at the lowest-order accuracy, it is conceivable 
that they may receive large modifications by including radiative and relativistic corrections. 
Moreover, for simplicity we have assumed the radiative correction to electromagnetic am- 
plitude does not introduce an imaginary part. One should realize this is just an (decent 



2 Besides the 1~0~ mode, other two-body decays of J/ip seem to also have a nearly orthogonal relative 
phase between a 7 and Qg , such as 0~0~ [15, 
in i/i' decays, the 1~0~ 4l| and _ 0~ mode 



;ay s ot J/tp seem to also nave a nearly orthogonal relative 
M J" 1 " [J HQ* 1+0 " and AW 0,0. Moreover 
42! . 43 1 seem also compatible with a large relative phase. 
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TABLE III: The Breit-Wigner, continuum and full cross sections (in units of fb) for e + e~ — > J/tj)+r} c 
at various T(nS) resonances. All the input parameters are the same as in Section [III Bl except a s 
is fixed to be 0.22. 



v 7 ^ (GeV) 


<7BW 


""cont 


Ofull 


9.460 


15678 


47.1 


14158 


10.023 


7165 


32.7 


6317 


10.355 


7948 


26.4 


7141 


10.579 


0.0026 


22.9 


22.5 



though) approximation. Despite this alertness, we still expect the qualitative feature, i.e., 
the large relative phase can withstand all these uncertainties. 



C. Continuum-resonance interference for double charmonium production 

For a given final state in e + e~ annihilation experiment near a vector meson resonance, it 
is always produced via two inseparable mechanisms- resonant decay and continuum produc- 
tion. A rough indicator about the relative strength of resonant electromagnetic amplitude 
to the continuum amplitude is characterized by 3B e + e -/a. For the first four T resonances, 
this factor is 10.2, 7.9, 8.9 and 0.0055 respectively. Therefore, for the three lower T res- 
onances, the J/ip + rj c production are dominated by the resonant decay, whereas for the 
T(4f>), which has a width about three orders of magnitude broader, one expects that the 
continuum contribution plays an overwhelmingly important role. 

We are interested to know the impact of the resonance-continuum interference on the 
observed cross sections. Assuming a 7 and a c differ by a Breit-Wigner propagator, one can 
express the full cross section near T peak as 



cm e e 



J/V> + Vc 



2 20 7T 2 eW|P 



X 



K 



9 s 
3a _1 \/s T 



s — M\ + *M T r tot 



(K + a 3g + 



(39) 



"991 1 



where T e + e - and r tot are the electric and total width of T. If the continuum term is dropped, 
this formula then reduces to the standard Breit-Wigner form: 



<7 B w e e 



12vrr e+e -r[T -> J/^ + rjc 

( " - My) 2 + m 2 r 2 



(40) 



'T 1 tot 



In Table II III we have enumerated various contributions to the J/ip + rj c cross sections at 
T(nS) peaks. One can clearly see the inclusion of the continuum contribution will reduce 
the peak cross sections by about 10% for the first three T states, whereas including the 
resonant contribution will reduce the continuum cross section by about 2% for T(45). This 
destructive interference can be attributed to the approximate 180° relative phase between 
a^g and a c . 
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10.52 10.54 10.56 10.58 10.6 10.62 



s (GeV) 

FIG. 5: The line shape of e + e~ — ► J/ifi + rj c near \fs = Myug\. 



The interference with continuum contribution also slightly distorts the Breit-Wigner 
shape of the production cross sections for the first three T resonances. However, one has to 
bear in mind that, for a thorough analysis, one has to carefully take the beam spread and 
radiative corrections into account, which requires lots of extra work and we leave them to 
the experimentalists. 

Thus far, the measured double charmonium production in B factories has been assumed 
to be entirely initiated by the continuum process, as represented in (152)1 . Experimentally, 
the resonant decay, despite its small magnitude, is encapsulated in the observed cross sec- 
tions. It is interesting to know how the line shape of J /if) + r] c near T(4S) peak would be 
affected by including this contribution. In Fig. (5) we have shown the various line shapes, 
with the contributions from several different sources juxtaposed. An interesting feature is 
that a dip is developed right on the T(4S') peak, which is again due to the destructive in- 
terference between the resonant strong decay and continuum amplitudes. Furthermore, we 
are reassured again that the radiative decay amplitude is unimportant. It will be great if 
someday experimentalists can do an energy scan and pin down this dip structure. To achieve 
this goal, the cross section must be measured very precisely, of course a very challenging 
task. We finally remark that, due to the aforementioned destructive interference, the true 
continuum cross sections should be slightly larger than the values quoted in (jMj) . which are 
in fact the full cross sections measured experimentally. 

IV. SUMMARY AND OUTLOOK 

In this work, we have performed a systematic study on T exclusive decays to vector 
plus pseudoscalar charmonium in NRQCD factorization framework. These exclusive decay 
modes can proceed via three-gluon, one-photon and two-gluon-one-photon, each of which 
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has been thoroughly analyzed. The relative phases among these amplitudes naturally arise 
as a consequence of the short- distance loop contribution. A particularly interesting finding is 
that the relative phase between strong and electromagnetic amplitude is nearly orthogonal, 
which is the same as that in various J /if} decay modes. 

The typical branching fractions of these decays are predicted to be of order 10~ 6 for the 
low-lying T(nS) states (n = 1,2,3). Future dedicated high-luminosity e + e~ facilities, e.g. 
Super B experiment, should be able to discover these decay channels readily. 

We have also investigated the impact of the continuum-resonance interference on the 
J /if> + T) c production cross sections at different T peaks. We find this interference will reduce 
the peak cross sections for the first three T states by about 10%. We predict there is a small 
dip in the line shape on the T(4S) peak. The current experiments are too rough to discern 
this delicate structure, perhaps the future Super B experiment can verify this prediction. 

A natural extension of this work is to investigate other exclusive double charmonium 
production processes from T decay. For example, T — > Xcj J /ip are particularly interesting 
channels to study, since the inclusive bounds for T — > Xcj + A have already been experimen- 
tally available |8j. Besides these double charmonium decay modes, one may also be tempted 
to apply the same formalism developed in this work to the processes T(J/ip) — > PV [44]. 
For the scarcity of theoretical investigations to these decay modes from the angle of pQCD, 
this study will offer us something worthwhile learning. Although it will no longer be as 
theoretically well grounded as the processes considered in this work, it should be viewed as 
an approach rooted in the time-tested constituent quark model, which has witnessed many 
phenomenological successes over years. 
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APPENDIX A: DERIVING ANALYTICAL EXPRESSION FOR / 

In this Appendix we illustrate how to simplify / effectively, so that we can obtain their 
analytical expressions. Repeatedly using kinematical relations stemming from the constraint 
k\ + k 2 = ®, plus fractional decomposition, we can reduce each fa in (1TTT) . (1T21 and (IT51) into 
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the sum of two-point tensor and three-point scalar integrals: 



h(0 



+ 



d 4 h f 2m 2 - 3m 2 c 

k x .[{2-0Q-P\ 



+ 



h ■ (3Q - P) 



mi 



^1^2 "a 1^2 '"'ht- 



ml. 



(k 2 - mD{k\ - h ■ P) k\{k\-k x -P) 



/ 2 (0 



d 4 h 



mi 



3m? 
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h.-P 



mi 
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(kl - m 2 )k 2 



+ 



h-ip-m 



mi 



(kl - m 2 )(k 2 - h ■ P) k 2 (k 2 -k 1 -P)_ 



MO 
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— 2m 2 c 
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PP 

r*2 
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(Al) 



(A2) 



(A3) 



(k 2 - m 2 )(k 2 - m 2 ) (kl-ml)(kl-ki-P) k\(k\ - k x ■ P) 



where £ = m 2 c /m 2 . While the two-point functions can be trivially handled, working out the 
three-point scalar integrals is more laborious but still straightforward. Here we just give 
their analytic forms: 



cm 



mi 



d 4 h 

'm 2 k 2 k 2 (k 2 -k x -P) 
2tanh _1 /3 ln£ + 2Li 2 



-2Li 2 



1 + P 



Li 2 
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(A4) 



(A5) 



It may be worth mentioning that if the well-known master formula for massive three-point 
scalar integral (i.e., equation (5.6) in (45[) is employed, one seems unable to obtain the 
correct expression for C 2 . To be specific, using that formula would render C 2 (|) = 0, which 
diametrically conflicts with the true value 4 In 2— 2ir/\^3. One can check our result is correct. 
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We now can express fa as follows: 
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Adding these three functions together then reproduces (TH|) and (|T5|) . 
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In this Appendix we illustrate how to reduce the one-loop four-point function in ( l20i) to 
the sum of simpler two- and three-point scalar integrals. With the aid of the kinematical 
identities arising from the constraint k\ + k 2 = P, we can disentangle this integral into three 
pieces: 
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Here we give the analytical expressions of two needed scalar 3-point integrals: 
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One then readily reproduces the analytic results shown in (12T1) and (|22[) . 
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